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ABSTRACT 

We report on the angular momentum content of heavily embedded protostars based on our analysis 
of the projec ted rotation veloci t ies (v sin is) of 38 Class I/flat spectrum young stellar objects recently 
presented bv lDoppmann et alJ ((2005) 1 . After correcting for projection effects, we find that infrared- 
selected Class I/flat spectrum objects rotate significantly more quickly (median equatorial rotation 
velocity ~ 38 km/sec) than Classical T Tauri stars (CTTSs; median equatorial rotation velocity ~ 18 
km/sec) in the p Ophiuchi and Taurus- Aurigae regions. Projected rotation velocity (v sin i) is weakly 
correlated with T e ff in our sample, but does not seem to correlate with Br 7 emission (a common 
accretion tracer), the amount of excess continuum veiling (r^), or the slope of the SED between the 
near and mid IR (a). The detected difference in rotation speeds between Class I/flat spectrum sources 
and CTTSs proves difficult to explain without some mechanism which transfers angular momentum 
out of the protostar between the two phases. Assuming Class I/flat spectrum sources possess physical 
characteristics (M»,R»,B») typical of pre-main sequence stars, fully disk locked Class I objects should 
have co-rotation radii with in their protostell ar disks that match well (within 30%) with the predicted 
magnetic coupling radii of I Shu et alJ l|1994fl . The factor of two difference in rotation rates between 
Class I/flat spectrum and CTTS sources, when interpreted in the context of disk locking models, also 
imply a factor of 5 or greater difference in mass accretion rates between the two phases. A lower 
limit of M ~ 10 -8 M Q /year for objects transitioning from the Class I/flat spectrum stage to CTTSs 
is required to account for the difference in rotation rates of the two classes by angular momentum 
extraction through a viscous disk via magnetic coupling. 

Subject headings: accretion, accretion disks — infraredrstars — stars:formation — stars:low-mass — 
stars:pre-main-sequence — stars:rotation 



1. INTRODUCTION 

Over the past two decades evidence has accumulated 
that the angular momentum content of low mass stars 
undergo significant modification during their first 100 
million years. This evolution has been best documented 
for cluster stars which have evolved to the T Tauri phase 
and beyond (for inf ormative reviews, see: iBodenheimerl 
119951 lMathieull200^1 . A key discovery is that the major- 
ity of low mass stars in the classical T Tauri phase are rel- 
atively slow rotators (v sin i < 20 km s _1 ; Bouvier et al. 
1986, Hartmann et al. 1986, Hartmann & Stauffer 1989, 
Bouvier 1990) contradicting expectations that accretion 
of high specific angular momentum material from a cir- 
cumste llar disk will drive stars to ro tate near breakup ve- 
locity ijHartmann &: Staufferlll989T) . Additionally, while 
the peak rotation velocity for a given population of young 
low mass stars consistently occurs at low velocities, a 
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high velocity tail of rapid rotators (v sin i > 50 km 
s _1 ) grows as a fraction of the total cluster population 
until an age of ~ 100 Myrs, based on o bservations of 
youn g clusters such as IC 2602 & 2391 (ISt auffer et all 
19971 iBarnes et al] 11999ft . a Per (IStauffer et all 119891 
Pressed 119921 I1994D and the Pleiades (IStauffer et al] 
19841 IStauffer fc Hartmandll987tlSoderblom et alJll993t 
Jones et alJll99ft iQueloz et al.lll998fl . 

The observational record is harder to interpret at the 
youngest ages. Investigations of Pre-Main Sequence 
(PMS) stars in Orion with photometric rotation periods 
have revealed a bimodal period distribution for higher 
mass (M > 0.25 M Q according to the models of D'Antona 
& Mazzitelli 1994, M > 0.5 M according to the models 
of Baraffe 1 998) stars, but only a si ngle peak for lower 
mass stars (Hcr bst et al1l2^Mj200l How ever, a simi- 
lar study of Orion PMS stars (|Rebullll2001^ was unable 
to detect this characteristic shift in the relation between 
rotation and mass. Likewise, early investigators found 
a link between rotation rates of T Tauri stars in Tau- 
Aur and Orion and infrared excesses (indicative of the 
presence of a circumstellar disk) such that slow rota- 
tors possesed large ex cesses, which rapid rotators lacked 
(jEdwards et al.lll993|) . These findings have been inter- 
preted in support of a paradigm in which T Tauri stars 
with large circumstellar disks (Classical T Tauri Stars 
[CTTSs], or Class II objects in the classification system 
of Lada 1987) are forced to rotate slowly via star-disk 
interactions, while T Tauri stars whose disks have dis- 
sipated (Weak-lined T Tauri Stars [WTTSs], or Class 
III objects) spin up to become rapid rotators as they 
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contract onto the main sequence. More recent investiga- 
tions in Orion, though, have found no such strong corre- 
lations between rotation rat es and photometric i nfrared 
excesses of T Tauri sta rs (jStassun et al.l 119991: iRebullI 
200I: lRhode et aUl200l . Indeed, it has even been sug- 
gested that the form of the correlation between rotation 
rates and infrared excesses that were reported in the past 
are contrary to what should be expected from theoreti- 
cal m odels of angular momentum transfer ijStassun et al.l 
120011) . Part of this confusion may stem from difficulties 
in using photometry alone to assess the presence or ab- 
se nce of disk material around T Tauri Stars; new work 
bv lSicilia-Aguilar et alJ l)2005l) using high resolution (R ~ 
35,000) spectroscopy to classify objects as either CTTSs 
or WTTSs on the basis of emission line widths suggests 
that photometric criteria for infrared excesses do not per- 
fectly segr egrate objects on the basis of their accretion 
strengths. iSicilia-Aguilar et al.l l)2005|) find that the use 
of spectroscopic criteria to separate CTTSs from WTTSs 
results in a clear difference between the rotation proper- 
ties of the two types of stars, where WTTSs appear to 
rotate significantly more rapidly than CTTSs. 

Even less is known about stars in earlier evolution- 
ary stages than T Tauri stars. To date, only a handful 
of embedded Class I and flat spectrum protostars have 
rotation rates measured fro m analysis of infrared pho- 
tosph eric spectral features QGreene fc Ladalll997i 1200(1 
12002]) . These objects appear to have faster rotation ve- 
locities (v sin i r~j 40 km s _1 ) than their more visible 
T Tau ri successors. In the optical. TWhite fc Hillenbrand! 
( 2004) (hereafter WH04) recently analyzed photospheric 
spectral features for a sample of 11 less heavily embed- 
ded (and thus optically visible) Class I sources in Taurus 
and found slower rotation velocities (v sin i ~ 20 km s _1 ) 
that were comparable to those detected in a companion 
sample of Class II Hcrbig-Haro drivers. It thus appears 
that Class I protostars can exhibit a considerable range 
of rotation velocities (10-50 km s -1 ) whose magnitude 
may be related to how deeply embedded the objects are. 
The lack of a single large, homogenous sample, however, 
has as yet prevented strong statistical constraints on the 
angular momentum distribution of Class I protostars. 

The observational difficulty (due to the high levels of 
extinction and significant veiling emission from circum- 
stellar material) of measuring rotation rates for Class I 
protostars is particularly unfortunate, as it appears a 
great deal of angular momentum evolution may occur 
during the protostellar accretion phase. The precursor 
angular momentum content of low mass stars can be as- 
sessed by observations of molecular cores within star for- 
mation regions, believed to be evolutionary predecessors 
to more evolved, collapsing protostars. Observations of 
these cores in molecular transitions which trace only the 
densest gas (e.g. NH 3 ; Goodman et al. 1993, Barranco & 
Goodman 1998) reveal linear velocity gradients expected 
for objects in solid body rotation, with a resultant esti- 
mated specific angular momentum of a typical dense core 
of 10 21 cm 2 s _1 . A comparison to the estimated specific 
ang ular momentum for isol ated T Tauri stars of 10 17 cm 2 
s _1 ijHartmann et alJl!986j) reveals that protostars must 
lose 99.99% or more of their angular momentum prior to 
entering the T Tauri stage. 

A number of theoretical mechanisms for transporting 
angular momentum from an accreting protostar have 



been suggested . Disk locking models l|Konig]| Il991t 
iShu et al.lll994f) envision a stellar magnetic field frozen 
into the midplane of a circumstellar disk at a point where 
the angular (Keplerian) rotation velocity is slightly lower 
than the stellar rotation velocity. This causes the star to 
be spun down by magnetic drag while the specific angular 
momentum of the disk material increases, later being dis- 
sipated by viscous transport or dis k winds. Wind models 
which do not invoke disk locking ijTout fc Pringkjll992ft 
involve a transfer of angular momentum from the star 
to circumstellar material through ionized stellar winds 
which are constrained to flow along rotating magnetic 
field lines. Other mechanisms, such as dynamical in- 
teract ions with gravitation al instabilities in protostellar 
disks (jDurisen et al.ll2~003j) or fragmentation of the core 
and disk into a multiple star system, also allow transfer 
of angular momentum away from the central star. 

Direct observational tests to distinguish between these 
angular momentum transport models have been difficult, 
particularly insofar as much of the evolution has already 
taken place by the time a protostar emerges from its ac- 
cretion envelope as a classical T Tauri star (CTTS). In 
an attempt to address this lack of observational data con- 
cerning the rotational properties of younger, more heav- 
ily accreting protostars, we present here an analysis of 
rotational velocities derived from a spectroscopic survey 
of Class I and flat spectrum protostars in the Tau-Aur, 
p Oph, Ser and o t her st ar forming regions presented by 
Dop pmann et al.1 l)2005f ). §2 describes briefly the obser- 
vations and projected rotation velocity (v sin i) deter- 
minations of the Class I and flat spectrum sample, as 
well as the creation of a comparison sample of T Tauri 
stars in Tau-Aur and p Oph with measured v sin i values. 
§3 explains the implementation of a method for numeri- 
cally inverting the observed v sin i distribution to reveal 
the most likely equatorial rotation velocity distribution. 
In §4 we compare the rotation of the Class I and flat 
spectrum sample to that of T Tauri stars, finding that 
the two samples are different to a high degree of statis- 
tical confidence. We analyze the implications of those 
differences in the context of various models of angular 
momentum transport for evolving protostars in §5, with 
§6 summarizing our main conclusions. 

2. OBSERVATIONS & PROJECTED ROTATION VELOCITY 
DETERMINATIONS 

2.1. Class I and flat spectrum sources 

The observed protostellar rotation velocities used in 
this study were derived from our large r Class I and flat 
spectr um survey described in full by iDoppmann et all 
(2005) (hereafter D05). We summarize here those de- 
tails of the observations and data analysis of the D05 
survey most relevant for this work. Spectra of Class I 
and flat spectrum (hereafter CI/FS) sources were ob- 
served over the course of 13 nights between May 2000 
and June 2003 u sing the near infrare d echelle spectro- 
graph NIRSPEC ijMcLean et alJl!998|) on the 10 m Keck 
2 telescope. Spectra were acquired with a 0.58" (4 pixel) 
wide slit (allowing R ~ 18,000) utilizing a 1024x1024 
pixel InSb detector array with the NIRSPEC-7 block- 
ing filter under median ~ 0.6" seeing. The NIRSPEC 
gratings were configured such that 4 prominent sets of 
spectral lines could be simultaneously observed in non- 
continuous orders: 1) the 2.1066 /im Mg and 2.1099 fim 
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Al lines, 2) the 2.1661 /xm Br 7 line, 3) the 2.2062,2.2090 
/itm Na doublet and the 2.2233 /im H2 line, and 4) the 
2.2935 ^m CO bandhead. 

All spectra were flat fielded, cosmic ray cleaned, sky 
subtracted, extracted, wavelength calibrated and cor- 
rected for telluric abs orption using standard IRAF pack - 
ages and procedures l|Massev et al.11199^: lMassevlll997j) . 
The observed NIRSPEC spectra of each protostar was 
then match ed to a grid o f synthetic spectra, constructed 
by MOOG <lSnedenll973H from the NEXTGEN model at- 
mospheres (Hauschildt et alJll999j) and modified to rep- 
resent various veilings, rotation velocities and radial ve- 
locities. Artificial rotation broadening of synthetic spec- 
tra was accomplished with the use of a rotation broaden- 
ing profile with an assumed limb-darkening coefficient of 
0.6 QOravll !)92 ;. Least squares minimization of the resid- 
uals between the observed spectra and synthetic spectral 
grid then produced estimates of the effective tempera- 
ture (T e ff), surface gravity (log g), radial velocity (v,.), 
projected rotation velocity (v sin i) and veiling (r^) for 
each protostar. The spectral resolution of the observa- 
tions translates to a velocity resolution of ~ 14 km s _1 ; 
detected rotational velocities were calculated with this 
slit broadening removed in quadrature. Error estimates 
for the derived rotational velocities were calculated from 
the shape of the local minimum in the fitting residual 
space sampled across the rotation velocity dimension and 
were typically 3-4 km s _1 . See D05 for individual v sin i 
velocities and uncertainties. 

To confirm the accuracy of our derived v sin i esti- 
mates, we have compared measurements of the rotation 
velocities of a number of our targets with values previ- 
ously reported in the literature for those objects. The 
derived v sin is for two spectral type standards in our 
sample (Gliese 791.2 and 1245a) were compared with the 
pr ojected rotation ve l ocities derived for the same objects 
bv lMohantv fc Basril |2003:) . We find rotation velocities 
for these stars of 34.3±1.3 and 22.5±2.2 km s" 1 respec- 
tiv ely, consistent within the errors with the values found 
bv lMohantv fc Basril p003h of 32±2 and 22.5±2 km s _1 . 
WH04 recently used optical spectra to produce v sin i 
estimates for 5 of the targets involved in this study: of 
those 5 shared targets, 4 (IRAS 04158+2805, DG Tau, 
GV Tau S and IRAS 04489+3042) produced derived ve- 
locities which matched within the errors on the measure- 
ments. WH04 also record an upper limit of 15 km s _1 for 
IRAS 04016+2610, while the measured v sin i value we 
find for IRAS 04016+2610 is significantly larger: v sin 
i = 46+2.7 km s _1 . WH04, however, flag their derived v 
sin i value for IRAS 04016+2610 as highly uncertain due 
to the extremely low signal to noise of their sprectra; as 
our spectrum is extremely well measured (S/N = 240) 
we assume the value reported by WH04 for this partic- 
ular object is in error due to the faintness of its optical 
detection. Our measured v sin i value of 24+2.0 km s _1 
for DG Tau (an optically visible T Tauri star with a flat 
SED) also matches well with the previously determined 
value of 21.7+6.3 km s _1 found bv lHartmann fc St auffcr 

Additionally, most field stars later than spectral type G 
have rotation velocities below the detection limit of this 
survey. As we would therefore expect, the 26 slowly ro- 
tating spectral type standards observed by D05 (all later 
than G7) provided only one with a significant (v > 14 km 



s _1 ) detection of rotation velocity: Gliese 569B. With a 
measured v sin i of 20 km s _1 , Gliese 569B has also been 
detecte d through adaptive o ptics imaging to be a binary 
system ijMartm et al J 120 00]) - it is possible the broader 
lincwidth of the composite spectrum is due more to com- 
bining spectra of individual components at different or- 
bital velocities than any intrinsic stellar rotation. The 
high level of agreeme n t betw e en the values repo r ted b y 
White fc Hillenbrand! |2004[) . iMohantv fc Basril l|2003ft . 
Hartmann fc Staufferl l)1989|) and our derived v sin i val- 
ues for shared targets, as well as our upper limits for 25 
of 26 late type field stars, support the accuracy of both 
our derived rotational velocities and our corresponding 
error estimates. 

The full CI/FS sample of D05 consists of 41 objects 
with measurable absorption lines. Some care, however, 
must be applied before blindly accepting all such ob- 
jects as input data for a study of this nature. Spectra 
of objects which are actually unresolved binaries could 
cause the characteristic rotational velocity derived here 
for protostars to be an overestimate, as the blending of 
spectral lines from two binary components at slightly dif- 
ferent radial velocities could be misinterpreted as rota- 
tional broadening of a single line in one object. To at- 
tempt to reduce this error, we have taken the following 
two steps: a) we have removed from our sample those 
objects known to be in binaries with separations on the 
order of or less than the slit width with which the ob- 
servations were taken (W L1 and Haro 6-28; Haisc h~et alJ 
20041 ISimon et al .11 995ft . and b) we have cross correlated 
our target spectra with similar spectral standards using 
the IRAF task FXCOR to see if any sources show struc- 
ture in their cross correlation results across multiple or- 
ders which suggests an observed spectra may be the sum 
of two individual components. This additional screening 
process identified only one object with signs of unresolved 
binarity; this source (IRS51) has been eliminated from 
our sample as well. An additional sign of binarity would 
be sources whose radial velocities are discrepant from the 
mean cluster velocity due to orbital motions. In our sepa- 
rate work inve stigating the radia l velocity distribution of 
these sources l|Covev et alJl200^ we have identified four 
objects within the CI/FS sample with radial velocities 
that are 3tr outliers with respect to the velocity of the 
local CO gas. Of these four sources, two are relatively 
slow rotators (v sin i < 30 km s _1 ) while two are more 
rapid rotators (v sin i > 40 km s" 1 ). We discuss more 
fully their possible binary status in our work on radial 
velocities, but with no concrete proof of their binarity, 
have left them in the rotation sample. Therefore, to the 
fullest extent possible we have tested for and eliminated 
unresolved binaries from our sample, but we welcome 
future adaptive optics observations of these systems to 
place tighter limits on the possible binary companions. 

Having been culled of possible binaries, the remain- 
ing 38 sources from D05 come from several nearby star 
formation regions, most notably Taurus Aurigae (10 tar- 
gets), p Ophiuchi (14 targets) and Serpens (9 targets). 5 
additional targets were selected from R CrA and Perseus. 
An additional 11 targets were observed in the course of 
the D05 survey but proved to have only smooth con- 
tinuum emission with no photospheric absorption lines 
in the spectral regions studied. All targets were cho- 
sen from surveys of local star formation regions utiliz- 
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FlG. 1. — Upper panel: A histogram, with 3 km/sec bins, of the 
observed v sin i values for a sample of 38 class I and flat spectrum 
protostars observed by D05. Lower panel: The distribution of spectral 
types for the protostars in the Class I and flat spectrum sample, based 
on the derived ^eff of the object and the dwarf T e f f -sp ectral type 
conversion presented by dc Jasicr & Nicuwcnhuiizcn (1987). 

ing IRAS, ISO, or ground based observations to assess 
the slope of membe r SEDs from the near to mid IR 
JWilking et allll989t IKenvon et all 1199(1 IGreene et al l 
1 1994 IKenvon fe Hartmannlll995l iBontemps et alJl200H 
iKaas et alJ I2004JT The sample was limited to objects 
with K < 11 in general, with some additional preference 
for objects with low bolometric luminosities to aggres- 
sively sample the low end of the mass spectrum. The 
top panel of Figure 1 displays a histogram of the raw 
v sin i velocities observed for the 38 protostars in the 
sample. 

To ensure the validity of conclusions concerning proto- 
stellar evolution derived from comparisons between the 
protostellar sample considered here and a sample of T 
Tauri stars (described in §2.2 below), it is important to 
ensure to the maximum extent possible that the two sam- 
ples contain objects with similar distributions of final 
main sequence stellar masses. Pre-main sequence stars 
largely contract along Hayashi tracks towards the low 
end of the main sequence, so the evolution of a low mass 
star from the CI/FS phase to the T Tauri phase should 
not introduce large changes in its effective temperature 
or spectral type in the absence of accretion. In order 
to aid the comparison of the rotation properties of our 
CI/FSs to those of a sample of T Tauri stars, we have 
calculated the approximate spectral type for each of the 
protostars in our sample based on the derived T e // of 
the object a nd the dwarf T P ff -spectral type co nversion 
presented bv lde Jaeer fc Nieuwenhuiizenl l)1987l) . In the 
range of T e //s relevant to this work, dwarf and PMS 
T e ff /spectral type relations differ only by a few hun- 
dred degrees K. Given the scale of other uncertainties in 
the process (see next paragraph) and that many T Tauri 



stars have been spectral typed using dwarf indices, we 
have chosen to perform our T e ///spectral type conver- 
sion with respect to a dwarf scale. The resultant ap- 
proximate spectral type distribution of the Class I and 
flat spectrum sample is presented in the bottom panel of 
Figure 1. 

CI/FS sources do, however, have derived mass ac- 
cretion rates M ~ 1CT 6 M yr" 1 ijWilking et all 119891: 
IGreene fc Ladal 12002): over an assumed lifetime of the 
CI/FS phase of 10 5 yrs, this would amount to an addi- 
tional maximum of 0.1 Mq being added to the source. 
Variations of 0.1 Mq in the mass of main sequence stars 
with spectral types typical of our CI/FS sample (~ K7) 
correspond to variations in spectral type of ~ 3 spectral 
type subclasses. Therefore, the unknown future accre- 
tion histories of CI/FS objects do introduce some un- 
certainty into our ability to use spectral type or T e ff 
as a proxy for the final main sequence mass of objects 
in our CI/FS sample. However, the scale of the result- 
ing bias does not seem large enough to create a signifi- 
cant mismatch between the final main sequence masses of 
CI/FS and T Tauri samples with identicial current spec- 
tral type ranges. We thus conclude that though more 
detailed knowledge of the future accretion history of our 
CI/FS sample would be desirable to properly construct 
an analogous sample of T Tauri stars, constructing sam- 
ples with similar current spectral type distributions allow 
a reasonable investigation of evolutionary effects. 

2.2. T Tauri stars 

In order to facilitate the interpretation of the measured 
v sin i distribution of Class I and flat spectrum sources 
in the context of Pre-Main Sequence (PMS) evolution, 
we have assembled a sample of ~ 150 pre-main sequence 
stars with spectr oscopically measured v sin i values from 
the lite rature llRebull et alJ | 2004| : |Doppmann et all 
20031 iWhite fc Basril l2Q03i I.Iavawardhana et alJ 



2003 IGreene fc Larlal 1200(1 ICla.rke fc BouvieH K 
Wichmann et al l 120001: iLuhman fc Riekel Il99£ 
a! 11992 LPrcibisch & Smith! 11992 [Padget 



Greene fc Ladalll997t iP reibisch & Smith 1997; Padgett 
1996; jKejiVjOji^llartmann 1995; Neuhaeuser et all 
" 19951 iGameiro fc Lagol 119931: iMagazzu et alJ I199S " 



w\)st iiiameiro & Laed litJHdt liviagazzu et alJ iwyA 
Bouyieilll99(l iHartmann fc Staufferl 119891: iWalter et al] 
1 9881 IBouvier et allll 9S(tI1 . As IClarke fc BouvieH iJSnOTO 



ler et all lU 

have found evidence for differing stellar rotation dis- 
tributions among T Tauri stars in Taurus and Orion, 
we have attempted to construct our sample exclusively 
from PMS in the same regions studied by D05. This 
restricts the PMS sample to Tau-Aur and p Oph since 
there are no well studied PMS stars in Ser and the 
D05 study includes only 5 objects in other regions, too 
few for statistical comparisons. We also restricted the 
comparison sample to stars of spectral type K0 and later 
to ensure that it sampled the same range of effective 
temperatures (T e ff) as the protostars in our sample 
(see Fig. 1). 

We have also divided these stars according to their 
pre-main sequence evolutionary phase; each star is clas- 
sified as either a Classical T Tauri Star (CTTS), a Weak 
T Tauri Star (WTTS), or Post T Tauri Star (PTTS), 
with classifications largely adopted from the primary ref- 
erence for an individual source. The notable exception 
concerns those objects detected as X-ray bright sources 
in the ROSAT All-Sky Survey <)Neuhaeuser et all 119951: 
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FlG. 2. — Display of the vsini and numerical spectral type distribution 
a given evolutionary phase, with the rotation velocity distribution for that 
up that sample shown on the right. 
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iWichmann eTaD [2000) ; given the uncertainty concern- 
ing the true ages o f these objects (jBriceno et all 119971 
iBriceno et al.lll999T) we classify them as PTTSs by de- 
fault unless they have been identified by other crite- 
ria (Ha equivalent width, IR excess, etc.) as a CTTS 
or WTTS. Projected rotation velocity distributions for 
these three comparison samples, along with the proto- 
stellar v sin i distribution shown in Figure 1, are shown 



(K0 = 0, K5 = 5, M0 = 10, M5 = 15). Each tier displays the sample for 
phase on the left, and the distribution of spectral types for stars making 



in Figure 2. Note that all samples (with the exception 
of PTTSs) contain objects with spectral types ranging 
from KO to mid M, and the median spectral type for each 
sample (again excluding PTTSs) is at the late K/early 
M boundary. While noting the accretion uncertainties 
raised in the previous section, we interpret the similar 
spectral type distributions of the various samples as a 
sign that each sample represents a relatively similar sam- 
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pling in final stellar mass. 

3. CALCULATING EQUATORIAL ROTATION VELOCITY 
DISTRIBUTIONS 

While analysis of spectroscopic lines can produce an es- 
timate of v sin i, without knowledge of the photometric 
period or other inclination diagnostics, it is impossible to 
solve the degeneracy between v and i for any particular 
star. Assuming that the rotation axes of a sample of stars 
with measured v sin i values are randomly oriented allows 
one to directly calculate the mean equatorial rotation ve- 
locity, v, from the mean of the observed values, y, where 
v = v sin i. As shown first bv iChandrasekhar fc Miinchl 
(1950), the means of the two distributions can be related 
as: 

4_ 

v = -y (1) 

7T 

IChand rasekha r~fe Miinchl l|1950(l also demonstrate that, 
in principle, the exact shape of the equatorial rotation 
velocity frequency function, f(v) (which we will refer to 
as the 'true' frequency function), can be derived by in- 
verting the relationship between f(v) and the observed 
rotation velocity frequency function, (j>(y). Accordingly, 
f(v) can be written as: 

K OV J v y2y/y2- V 2 

which can be numerically calculated, as shown bv lGaigel 
given a continuous form of cj>(y). The remainder 
of this section merely describes our implemen tation of the 
method originally presented by Gaige (1993): interested 
readers are encouraged to consult the primary work for 
a fuller explanation. 

Our observations are simply the rotation velocities of 
a discrete sample of stars. These data are inherently 
dis-continuous, even when formed into a histogram as in 
Figure 1. To translate this set of N discrete observations 
into a continuous frequency function, <fr(y), we sum a set 
of gaussians centered on each observed rotation velocity 
(Yi) with a width set by the uncertainties in the derived 
velocities (er = 4 km s -1 ), as shown in Equation 3: 

1 1 9 9 

^) = ]vE^(^+^) (3) 

With a continuous observed frequency function in hand, 
we are now able t o apply the n umerical inversion algo- 
rithm presented bv iGaieel J1993) in their Equation 35 to 
solve our Equation 2 and calculate the 'true' frequency 
function, f(v). 

To test that we have accurately implemented this nu- 
merical inversion technique, we have created simulated 
observed velocity distributions as inputs for this inver- 
sion method, and compared the results with the underly- 
ing 'true' equatorial velocity distributions originally used 
in generating the simulated observed sample. To create 
an simulated observed velocity distribution as input for 
the numerical inversion algorithm, we first randomly se- 
lected a sample of N stars from a given equatorial rota- 
tion velocity distribution. Each star was then assigned a 
random inclination, selected in a manner consistent with 
an isotropic distribution of rotational axes, and the resul- 
tant simulated observed v sin i for that star is calculated. 




FlG. 3. — A test of the accuracy of the numerical inversion tech- 
nique used to determine the 'true' distribution of equatorial rotational 
velocities from a set of observed v sin i values. These tests were run on 
a sample of 1000 synthetic stars with velocities distributed uniformly 
between and 60 km s~ 1 and random inclinations. Top panel: A 
comparison between a histogram of the simulated 'observed' vsini val- 
ues and the smoothed, continuous rotational frequency function, 4>(y). 
Middle panel: A comparison between a histogram of the input equa- 
torial rotational velocities and the numerically inverted determination 
of the 'true' distribution function. Bottom panel: a comparison of the 
cumulative frequency functions for both the input equatorial velocity 
distribution (solid) and the derived 'true' cumulative frequency func- 
tion (dashed). 

All N artifically observed velocities were then combined 
into an aggregate velocity distribution; this distribution 
was analagous to the distribution of velocities from D05 
which form the raw observational data for this study. 

To complete our test of the inversion method, we then 
completed the necessary steps to derive a 'true' equato- 
rial rotation velocity distribution from an observed set 
of data: a continuous frequency function was formed 
from the discrete set of simulated observed velocities, 
and this continuous frequency functio n was then nu meri- 
cally inverted following Equation 35 of Gaige ( 1993]). We 
tested the method with two types of input velocity dis- 
tributions: a sample of 1000 stars with a single allowed 
equatorial velocity (v eq = 60 km s _1 ) and a sample of 
1000 stars where equatorial velocities were randomly dis- 
tributed with a uniform probability within the range 
km s _1 < v < 60 km s . The basic form of the in- 
verted 'true' distribution is correct in each case, deriving 
a single velocity population for the first case and a more 
broadly distributed population in the second case (shown 
in Figure 3). All differences between input and output 
distributions were consistent with the gaussian smooth- 
ing and finite statistical noise of the input samples. 

The middle panel of Figure 3 shows the amplitude of 
statistical fluctuations in both the input equatorial rota- 
tion velocity distribution and the inverted 'true' distri- 
bution. The input distribution was selected such that all 
velocities had a uniform probability, so deviations from 
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a constant number of objects per bin represent statis- 
tical fluctuations in the sampling of the underlying dis- 
tribution, consistent with the level expected by Poisson 
statistics. This random velocity sampling, coupled with 
each star's randomly assigned inclination, results in some 
level of disagreement between the 'observed' sample and 
an ideally sampled, perfectly isotropic distribution of ro- 
tational axes required for the method to return the cor- 
rect result. Inspection of the middle panel of Figure 3 
indicates that, even for well sampled populations, devi- 
ations between the inverted 'true' distribution and the 
probability function from which stars were selected may 
be as high as ±15% at any individual velocity. 

The cumulative velocity distributions for both the in- 
put and inverted 'true' distributions are shown in the 
bottom panel of Figure 3. These distributions, which 
display the integral of the probability distribution from 
to a given velocity, appear much less noisy to the eye. 
Individual peaks and valleys from the middle panels are 
smoothed over by the integration, resulting in a good de- 
gree of agreement between the input and inverted 'true' 
distributions, within the limits of the resolution of the 
technique. 

Having tested the accuracy of the inversion routine in 
the well sampled regime, we also simulated the results ob- 
tained from sparsely sampled simulated observed distri- 
butions, to test the statistical fluctuations which could be 
introduced into our method as a result of using an input 
sample of only 38 stars. This test shows the same types 
of errors as in the N= 1000 test, but amplified, as statis- 
tical fluctuations in the input sample are correspondingly 
larger as a fraction of the number of stars within a given 
range of velocities. The underlying form of the derived 
'true' distribution still appears fundamentally correct, as 
a broadly distributed population, but statistical fluctua- 
tions in the resultant probability distribution appear to 
have increased to ± 30% at any given velocity. However, 
the cumulative frequency function continues to integrate 
over many of the sharp peaks and valleys in the inverted 
probability distributions, and the agreement between the 
two functions appears good, even with only N=38 objects 
to base the inversion on. 

We conclude from these tests that the underlying in- 
version method is sound and has been implemented cor- 
rectly, though our sparsely sampled input distribution 
will clearly affect the accuracy of our derived 'true' dis- 
tribution. Thus, in Figure 4 we present the numerical in- 
version of the observed v sin i distribution derived from 
the 41 protostars in our sample. We expect from our 
previous tests that statistical fluctuations on the order of 
±30% at any given equatorial rotation velocity will exist 
in our derived probability function (middle panel). Our 
cumulative frequency function (lower panel), however, in- 
tegrates over many of these fluctuations and presents a 
more statistically robust assesment of the distribution of 
equatorial rotational velocities for CI/FS protostars. 

For comparison, Figure 5 presents the results of this 
numerical inversion technique as applied to our sample of 
CTTS rotational velocities obtained from the literature. 

4. ANALYSIS 

4.1. Statistical tests of angular momentum evolution 
and possible sample biases 




FlG. 4. — The inversion of the v sin i distribution measured by D05. 
Top panel: a comparison between the observed v sin i distribution 
(dashed histogram) and the resultant continuous frequency function, 
4>{y). Middle panel: The inverted 'true' frequency function, show- 
ing the probability of a class I or flat spectrum protostar having a 
given equatorial rotation velocity v. Bottom panel: The cumulative 
frequency function of the 'true' equatorial rotation velocity distribu- 
tion of class I and flat spectrum protostars. 




FlG. 5. — An inversion of the observed fsim frequency function to 
produce the 'true' frequency function, as in Figure 4, but for the CTTS 
comparison sample. 
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A first test for the presence of angular momentum evo- 
lution between the CI/FS sample and that of more vis- 
ible, and presumably more evolved, PMS stars can be 
determined from the observed v sin i distributions con- 
tained in Figure 2. The samples of observed v sin is 
taken from the literature for each T Tauri stage can be 
compared with the observed v sin is of the CI/FS sample 
via a two sided Kolmogorov-Smirnov (KS) test for sim- 
ilarity. The KS test compares the cumulative frequency 
functions of the observed rotation velocity distributions 
for the two samples; the maximum distance between the 
two cumulative frequency functions can be used to cal- 
culate the probability th e distributions we re drawn from 
the same parent sample l|Press et~aTllT986l) . The KS test 
reports a 1.4x 1(T 10 , 1.6x 10" 5 and 1.6 x 10" 4 probability 
that the CTTS, WTTS and PTTS samples, respectively, 
were selected from the same parent distribution as the 
CI/FS sources. 

The KS tests simply confirm statistically what can 
be seen by comparing the distributions in Figure 2 by 
eye: the mean rotation velocity of the CI/FS sample is 
considerably higher than that of the PMS samples, and 
the shapes of the two distributions differ as well. The 
PTTSs are the most similar of the PMS samples to the 
CI/FS sample, a fact attributable to the higher rotation 
rates reached by the PTTSs as they contract onto the 
main sequence with little angular momentum loss. To 
investigate the amount of angular momentum evolution 
which must have taken place between the CI/FS and the 
CTTS phases, we select as representative velocities of 
each phase the velocity at which the 'true' cumulative 
velocity distribution reaches the 50% level. These veloc- 
ities are ~ 38 km s _1 and ~ 18 km s _1 , respectively. 

To seek evidence for the nature of any protostcllar an- 
gular momentum regulation mechanism at work during 
the CI/FS phase, we have searched for correlations be- 
tween projected rotation velocity and other observable 
protostellar parameters indicative of their evolutionary 
state. In doing so, we have found no dependence of pro- 
jected rotation velocity on accretion rate as measured by 
either Br 7 emission or excess continuum emission (r/-; 
see Figure 23 by D05). There also appears to be no de- 
pendence of rotation velocity on the slope of the SED 
(a) between 1 and 10 /jm for Class I or fiat spectrum 
sources. 

The only correlation which does appear to be present 
is one between projected rotation velocity and T e ^/, as 
shown in Figure 6. The correlation has a large amount of 
scatter, represented by the fairly low (0.53) correlation 
coefficient. Removing the two high temperature, high v 
sin i objects further weakens the correlation (correlation 
coefficient = 0.33), though the slope of the best fit line 
appears relatively similar between the two fits. 

It is important to note that both the CI/FS and CTTS 
samples suffer from observational biases which could 
skew the resultant v sin i distributions. In the case of 
the CI/FS sample, binarity (tested for and previously de- 
scribed in section 2.1) and veiling biases are likely present 
at some level within the sample. A veiling bias will re- 
sult from difficulty in detecting photospheric absorption 
lines from CI/FS sources whose spectra is affected by 
strong rotational broadening or excess continuum emis- 
sion at 2/xm from circumstellar material. Both the broad- 
ening due to rapid rotation and the veiling caused by ex- 
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FlG. 6. — Top: Teff vs. vsini for Class I/flat spectrum sources 
arc plotted as diamonds for individual sources. The best fit linear 
least squares relationship is ovcrplotted as a line; the raw correlation 
coefficient of the data, as well as the coefficients of the best fit line arc 
displayed within the plot. Bottom: As above, but after removing the 
two high temperature, high vsini objects from the sample. 

cess continuum emission make photospheric lines appear 
shallow, and thus harder to detect in a spectra with a 
finite signal-to-noise ratio. Therefore, the fact that 11 
of the 52 Class I/flat spectrum sources targeted for the 
survey of D05 did not reveal detectable photospheric ab- 
sorption lines indicates that those sources were either a) 
rapid rotators, b) heavily accreting sources, or c) ear- 
lier than G type, where the photospheric absorption fea- 
tures analyzed here become too weak to detect. The last 
possibility, that these stars are merely hotter than those 
which show photospheric absorption lines at 2 /xm, can 
be ruled out for most of the sample on the basis of a 
lack of a strong HI Br 7 absorption line, a common spec- 
tral feature of hot stars. None of the stars for which we 
fail to detect photospheric absorption features possessed 
a Br 7 absorption line, and 9 of the 11 had Br 7 emis- 
sion features in excess of 1 A equivalent width, widely 
interpreted as a sign of active accretion. As a result, 
the detections of D05 are incomplete at the 20% level 
and likely biased against the most heavily accreting or, 
most significantly for this work, rapidly rotating proto- 
stars; see Figure 23 of D05 for the region of yk and v 
sin i parameter space to which the survey is sensitive. 
Our characteristic velocity for this phase is thus likely 
an underestimate of the median velocity of a truly com- 
plete sample of CI/FS sources in the absence of binary 
contamination. 

The CTTS sample is biased in the opposite sense, such 
that our estimate of the characteristic velocity of this 
phase is an overestimate. This bias results from the ob- 
servational limits of the surveys used to measure many 
of the stars within the sample, which were only sensitive 
to rotational velocities > 10 km s _1 . Thus, the large 
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peak in the CTTS distribution at 10 km s _1 is formed 
with a number of stars whose v sin i detection is merely 
an upper limit; replacing those velocities with the true, 
and presumably lower, rotation velocity would lower the 
characteristic velocity we have determined for the CTTS 
phase. 

As the CTTS sample is biased towards rapid rotators, 
and the Class I/flat spectrum sample is biased towards 
slower rotators, the ~ 20 km s _1 difference between the 
characteristic equatorial velocities of each phase repre- 
sents a lower limit to the true change in equatorial rota- 
tion velocity evolution between the two phases. 

4.2. Comparison with previous work 

In a recent study of optical spectra of 11 Class I sources 
and 21 Class II sources (all known to be driving Herbig 
Haro outflows) in Taurus, WH04 found no evidence for 
differences in the rotation properties of the two samples. 
As IR selection of Class II objects and optical selection 
of CTTSs are thought to identify identical populations 
of pre-main sequence stars, WH04 appear to draw an 
opposite conclusion to that found in §4.1. Deriving ro- 
tation velocities from analysis of optical spectra show- 
ing photospheric absorption features, WH04 find nearly 
identical means for the two samples. The mean of the 
Class I sample was < v sin i > — 18.95 km s _1 while 
the mean of the Class II sample was < v sin i > = 17.39 
km s _1 such that a KS test revealed that there is 7.9% 
of both samples sharing the same parent distribution. 
The Class II sample does return a lower median v sin 
i (14.4 km s^ 1 ) than the Class I sample (18 km s _1 ). 
As shown previously in §2.1, v sin i determinations for 
the 5 stars in common between this program and that of 
WH04 agree well: 4 in 5 have determinations identical 
within the errors and another object (IRAS 04106+2610) 
is only marginally detected (S/N = 2.1) by WH04 with 
a v sin i value flagged as being highly uncertain. Thus, 
both samples appear to agree well on v sin i for common 
objects with observations in the high S/N regime. 

While measurements of v sin i for individual sources 
do appear to match well for those objects observed by 
both programs, the mean projected rotation velocity for 
Class I objects in Taurus derived by each program does 
differ (WH04: 19.1 km s" 1 , this work: 30.1 km s" 1 ). 
This difference appears to be due to an inverse correla- 
tion between the rotation rate of the Class I object and 
its optical visibility. Three Taurus Class I objects for 
which we derive high rotation rates from well detected 
infrared absorption features (IRAS 04016+2610 [46 km 
s" 1 ], IRAS 04295+2251 [51 km s" 1 ] and L1551+IRS5 
[31 km s -1 ]) appear to be so faint in the optical that 
WH04 were unable to derive confident estimates of v sin 
i. Accordingly, we conclude that the difference between 
the Class I v sin i distribution presented here and that 
of WH04 can be attributed to the fact that the fastest 
rotating Class I sources appear to be optically dim, likely 
due to higher levels of extinction, consistent with our pre- 
vious finding that more heavily embedded Class I sources 
have higher mean v sin i values than more optically vis- 
ible Class II/ CTTS sources 5 . 

5 We note that protostellar envelopes may be flattened along 
their rotation axes, either due to rotation itself or from clearings 
produced by bipolar jets. If this is the case, extinction will be 



4.3. Rotation velocity comparisons between regions 

The bulk of the information currently known con- 
cerning the rotation velocities of CTTSs has resulted 
primaril y from studie s of two star formation re gions : 
Taurus feouvier et alJfl986t lHartmann fc Stauffeilll989l: 
Bouy erl I1990D and Orion iStassun et all 119991: iRebullI 
200lURhode et al]l200lt ISicilia-Aguilar et alJl2005l) . Re- 
cent work bv lClarke fc Bouvierl l)200ff ) suggests that ro- 
tation properties for samples of stars selected as CTTSs 
in each region may be significantly different, while 
ISicilia-Aguilar et alJ l|2005ft find a typical CTTS rotation 
velocity in Orion which agrees well with our comparison 
sample of CTTSs drawn from Taurus and p Ophiuchi, 
implying a relatively small role for environmental ef- 
fects on the angular momentum evolution of young stars. 
With measurements of rotation velocity for ~ 10 Class 
I protostars in each of three star formation regions (p 
Ophiuchi, Serpens and Taurus), it is possible to begin 
to probe the degree to which the rotational properties of 
Class I sources differ across regions. 

Of the three regions, Serpens has the largest mean ob- 
served v sin i in our sample (36.8 km s _1 ), followed by 
p Ophiuchi (31.1 km s -1 ) and finally Taurus (30.1 km 
s^ 1 ). Two sided KS tests indicate Taurus and Serpens 
differ most strongly, but still have a 18% chance of shar- 
ing a parent distribution, so these differences between 
rotational properties of Class I objects in differing star 
formation regions are not yet statistically significant. 

5. ANGULAR MOMENTUM DISCUSSION 

The total angular momentum (L tot ) of a star can be 
expressed as L to t — where I* is the stellar mo- 

ment of inertia and u> is the angular velocity of the star. 
Assuming spherically symmetric, solid body rotation of 
a protostar with a conn ective interior, I* ~ 0.2M*Rl 
ijClaret fc Gimenedll990j) . and u = Thus, the an- 
gular momentum content of a given star can be simply 
written as Ltot ~ 0.2M*R*v. This functional form pro- 
vides a simple physical framework in which to interpret 
the implications of the difference in the median equato- 
rial velocities for both the CTTS and the CI/FS samples 
under various theoretical frameworks for protostellar evo- 
lution. 

5.1. Evolution With No Extraction of Protostellar 
Angular Momentum 

For there to be no extraction of angular momentum 
from the protostar between the CI/FS and CTTS phases, 
the difference of a factor of two in the characteristic me- 
dian equatorial velocities requires that the product of 
M*R* must be double that of the CI/FS stage by the 
CTTS phase. Absent any significant changes in the inter- 
nal structure of protostars between the two phases (such 
as the onset of dueterium shell burning, expected to in- 
flate the stellar radius by a factor of two or greater but 
only occur once stars have grown to > 3 M Q ; Palla 

1999), R* grows slowly as the total mass of the star in- 
i 

creases: R* cx M* . Stellar radius is believed to decrease, 
not increase, as PMS stars evolve toward the main se- 
quence, but this is not a significant effect in the evolution 

inclination dependent, and a bias towards optically visible sources 
will result in a bias towards sources with small sin is, and thus 
small v sin is. 
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from the CI/FS phases; D05 find that the CI/FS sample 
had surface gravities similar to those of CTTSs. Thus, 
an increase in M*R* by a factor of two and no signifi- 
cant changes in internal structure implies an increase in 
M* by a factor of ~ 1.7. Assuming a characteristic mass 
of 1 Mm for our sa mple, as predicted by the models of 
iBaraffe et aD l(l99^ for a T eff = 3800 K star at 1 Myr, 
an increase of this size implies accretion of 0.7 M Q of 
material. 

CI/FS objects are commonly thought to be nearing the 
end of the main mass accretion phase (a phase fully com- 
pleted by the Class II/CTTS stage), which should rule 
out a near doubling of the protostellar mass. This con- 
clusion is butressed by the comparison of spectral types 
for the CI/FS and CTTS samples (as seen in Figure 2 
and discussed in §2), consistent with little to no signifi- 
cant additional mass accretion between the two phases. 
However, even if mass accretion takes place at a signifi- 
cant level between the CI/FS phase and the CTTS phase, 
allowing that additional mass to lead to a decrease in 
rotation velocity requires that the newly accreted mate- 
rial have significantly lower specific angular momentum 
than that whic h was previously accreted by the proto- 
star. iPophaml l)1996j) suggest that FU Ori-like accre- 
tion events can result in mass accretion distinguished 
by markedly lower specific angular momentum than that 
which is accreted during more quiescent periods of proto- 
stellar evolution, precisely the effect needed to drive the 
rotation velocity of CTTSs down from the CI/FS level. 
However, the amount of FU Ori accretion required would 
be unusual given current estimates of the FU Ori duty 
cycle: pha ses of rapid mass accretion at the 10~ 4 M Q 
yr -1 level l)Hartmann fc Kenvonlll99fil) would require a 
cumulative total of ~10 4 years to be spent in the high 
accretion state to increase the protostellar mass by ~1 
Mq , a factor of 20 greater than the current estimated FU 
Ori duty cycle of ~500 years in the high accretion state 
for each low mass star. 

Ultimately, a model where the change in equatorial ro- 
tation velocity observed between the Class I and CTTS 
phases is due solely to redistribution of the current proto- 
stellar angular momentum is difficult to construct. This 
model is also difficult to reconcile simultaneously with 
current measurements of both Class I and FU Ori accre- 
tion rates and lifetimes, casting doubt on the assumption 
that no extraction of angular momentum takes place dur- 
ing the Class I phase. 

5.2. Deriving J assuming efficient extraction of angular 
momentum 

The opposite limiting case to that presented in §5.1 
would hold that the difference in v sin is of the two sam- 
ples is due completely to extraction of angular momen- 
tum from the protostar. We thus investigate the amount 
and rate of angular momentum extraction implied under 
a model where the physical parameters (M„, R») of a 
typical CI/FS object identically match those of a typical 
CTTS. Assuming canonical values of M* = 1.0 M Q and 
R* = 2.7 R© and expressing the change in angular mo- 
mentum between the two phases as AJ= 0.2M»i?*Ai>, 
we find AJ = 1.5 x 10 50 gm cm 2 s _1 . Further assum- 
ing a mean age of 10 5 years for our CI/FS sample and 
10 6 years for our CTTS sample implies ~ 0.9 Myrs of 
elapsed time between the Class I and CTTS phases, and 



allows a calculation of J = = 5.3 x 10 36 gm cm 2 s 2 . 

This derived value for J assumes no particular model for 
how angular momentum extraction is performed, but de- 
pends only on a) the assumed separation in ages between 
the Class I and CTTS phases [.9 Myrs], b) the assumed 
values for M„ and R* [M, = 1.0 M , R* = 2.7 R ], 
as well as assuming those values characterize both the 
CTTS and CI/FS phases, and c) the derived difference 
[~ 20 km s" 1 ] between typical CI/FS and CTTS equa- 
torial rotational velocites. 

5.3. Extraction of angular momentum via disk locking 

Disk locking models describe extraction of protostellar 
angular momentum via magnetic coupling of the proto- 
star through accretion columns to a region in the disk 
whose Keplerian angular velocity is slightly lower than 
the angular velocity of stellar rotation. Drag on the stel- 
lar magnetic field lines as they rotate through the disk 
then begins to slow the stellar rotation until the star, and 
likewise its magnetic field, are rotating with the same an- 
gular velocity as the disk. Once the protostar and cou- 
pling radius of the disk are rotating with the same an- 
gular velocity, the magnetic fields are stationary in the 
rotating reference frame of the disk - no magnetic drag 
occurs and no additi onal angular momentum is ext racted 
from the protostar l|Shu et al.lll994t lKonig]lll99lD . Ob- 
servational evidence of excess blue continuum emission 
(assumed to arise from accretion shocks; Bertout, Basri 
& Bouvier 1988) and the presence of magnetic fields 
strong enough to co uple protostars to the i onized por- 
tions of their disks l)Johns-Krull et aTl 1 1999(1 have been 
interpre ted in support of this p aradigm. However, recent 
models l)Matt & Pu dritz 2004) which attempt to account 
for dynamical evolution of magnetic field topologies in 
response to the star disk interaction appear to have diffi- 
culty extracting angular momentum at the rates assumed 
necessary to account for the slow rotation rates (~ 10% 
of breakup) observed in CTTSs. Additionally, observa- 
tional evidence f or the correlations predicted by these 
theories is mixed (|.Tohns-Krull fc Gaffordll2002T) . Despite 
these difficulties, disk locking theories present the most 
fully developed theoretical framework within which to 
explore the implications of the rotation velocity evolu- 
tion of protostars as they pass from the CI/FS into the 
CTTS phase. 

5.3.1. The typical Class I co-rotation radius 

We determine here the corotation radius of the pro- 
tostellar disk to which a typical CI/FS object would be 
magnetically coupled if fully disk locked. We note that 
the corotation radius will be a lower limit to the radius at 
which the protostar couples to the disk, as the time scale 
over which disk locking removes angular momentum from 
a protostar may be sufficiently long that these proto stars 
are still being magnetically braked (Hartmann 2003), re ~ 
quiring coupling radii outside the co-rotation radius. It 
should be noted, however, that differential rotation be- 
tween the star and disk should drive significant changes 
in the magnetic field topology, leading to reconfigurations 
that lessen the efficienc y of angular momentum transport 
UMatt fc Pudrit^TOl. 

To derive the typical co-rotation radius implied for the 
38 km s~ 4 median velocity of CI/FS sources, we sim- 
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ply invoke Kepler's third law, assuming circular disk or- 
bits where the mass of the disk has a negligable effect 
on the gravitational potential of the system, and thus 

Plisk = (fir 7 " 3 - ^ s we see k ^ e ra dhis of the disk where 
the Keplerian period is equivalent to the stellar rota- 
tion period, we can thus express the corotation radius 

as r co = ( — 1 ) 1 ^ 3 - Assuming a canonical mass of 
0.5 M© and radius of 3.1 R© for our CI/FS sample (we 
adopt 0.5 M© and 3.1 R© here instead of the 1.0 M© 
and 2.7 R© assumed in previo us sections to facilitat e 
comparison with the results of iGreerie fc Ladal (1200^) : 
assuming a mass of 1.0 M© and radius of 2.7 R© here 
would lead to an increase in the derived radius by 15%) 
implies a co- rotation radius of 8.6 R©, or 2.8 R*. We 
note that this value is close to the 2.1 R* magnetic cou- 
pli ng radius derive d for the Cla ss I object YLW 15 A 
by IGreene & Ladal i200l . The IGreene fc Ladal ll200l 
result was determined by simultaneously calculating the 
mass accretion rate and magnetic coupling radius nec- 
essary to account for the observed acc retion lum i nosity 
using the magnetic accretion models of lShu et alJ l)1994f) 
with protostellar parameters identical to those presented 
above. The agreement between these parameters sug- 
gest that the location at which co-rotation occurs in a 
typical CI/FS object matches well with the location at 
which coupling would lead to the observed rotation rate 
of CI/FS objects if magnetic torques are strong enough 
to lock the protostellar disk to the protostar. 

5.3.2. Implied protostellar mass accretion rates 

Magnetospheric accretion models typically assume 
that the circumstellar disk is truncated at the Alfvcn 
radius, where the energy density of the stellar magnetic 
field is matched by the kinetic energy density of the in- 
falling material. This truncation radius (r t ) can be ex- 
pressed as 



r t — R*( 



BtRl 



(4) 



'GMM 2 

It is also near this truncation radius that magneto- 
spheric accretion models typically locate the co-rotation 
point to which the star is locked, such that = tJdisk- 
Thus, for a star which is fully disk locked to the Keple- 
rian orbital velocity at the disk truncation point, we can 
express the stellar rotation velocity as 



v* = vdiskin) 



R* 
n 



IGMRI 



(5) 



Combining equations 4 and 5 allow us to express v* in 
terms of fundamental protostellar parameters, such that 



G 5 M$M 3 , 



(6) 



Therefore, viewed in the context of magnetospheric ac- 
cretion, along with the assumption that the basic stellar 
parameters (M*, R*, B») do not change significantly be- 
tween the CI/FS and CTTS phases, the observed change 
in the mean rotation velocity of the two samples would 
naturally be attributed to a change in the mean pro- 
tostellar accretion rates between the two phases if the 
protostars are fully disk locked. As v» cx M?, a decrease 



by a factor of 2 in the characteristic rotational velocities 
from the CI/FS to CTTS phase would represent a factor 
of 5 decline in the typical accretion rate from the CI/FS 
to the CTTS phase. This finding fits well into the ob- 
servational paradigm which assumes CI/FS objects are 
younger, more heavily accreting protostars than CTTSs. 
We must note, however, that the stellar rotation velocity 
is more weakly dependent on M than any other stellar 
parameter: thus, the degree to which stellar parameters 
change between the CI/FS and CTTS phases must be in- 
vestigated before protostellar rotation velocity can serve 
to test the validity of magnetospheric accretion models 
by comparing rotation velocities to observed accretion 

rates. 

\ noted by[Hartmann ( 2002), models which assume a 
transport of angular momentum away from a protostar 
through a viscous circumstellar disk require inward ac- 
cretion of mass in the disk to allow the outward transport 
of angular momentum. We can thus attempt a crude es- 
timation of the typical mass accretion rate between the 
CI/FS and CTTS phases by estimating the mass accre- 
tion rate necessary to allow the J derived in §5.2. We 
equate the rate of change of the angular momentum of 
the protostar with the rate at which angular momentum 
is transported away from the co-rotation radius via accre- 
tion through the disk. Following Hartmann's equation 5, 
but ignoring the spin-up torque of the additional newly 
accreted material to derive a lower limit to the needed 
accretion rate, we find: 



J « Muj(r co y co = M-^r 2 co 
-ft* 



(7) 



Applying a corotation radius of 9.9 R© derived in §5.3.1 
for a 1 M© protostar with a 2.7 R© stellar radius (R*) 
and a 28 km s _1 stellar velocity chosen to be intermediate 
to the 38 km s~ x mean CI/FS velocity and the 18 km s _1 
mean CTTS velocity, we can solve for the M applicable 
to objects transitioning from the CI/FS phase to CTTSs. 
We note that this calculation implicitly assumes that M 
decreases slowly enough that the star stays magnetically 
coupled to its disk over the CI/FS to CTTS transition pe- 
riod, allowing continuous transfer of angular momentum. 
We therefore find a lower limit to the mass accretion rate 
between the CI/FS and CTTS phases of M - 10~ 8 M© 
yr _1 would be required for a viscous disk to transport 
the angular momentum our observations require to be 
extracted between the two protostellar phases under the 
assumptions of §5.2. This lower limit fits well with the 
accretion rates measu red for CTTSs (~ 10 -8 M© yr" 1 ; 
IMuzerolle et alJ|1998|) . This accretion rate lower limit 
would increase if the age gap between the CI/FS and 
CTTS phases were lowered from the assumed value of 
0.9 Myrs, or if the effects of the spin-up torque due to 
accreted material were accounted for; either of these ef- 
fects would drive the intermediate accretion rate closer 
to the 10 -6 M© yr _1 estimated for the protostar YLW 
15 A in the context of magnetospheric accretion models 
JGreene fc LadalEfiol. 

Lastly, iHartm ann ( 2002J) also notes that the timescale 
over which the disk-locking process removes angular mo- 
mentum from a young star is a dynamic quantity, sensi- 
tive to changes in the star's rotation velocity, magnetic 
field strength, mass, and mass accretion rate. Specifi- 
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cally, iHartmaiml l)2002() predicts in his equation 8 a disk- 
locking timescale which varies as: 

tdb > 4.5 x 10 6 yrM . 5 MZ^f, (8) 

where Mo. 5 js the stellar mass in units of half a so- 
lar mass, M_8 is the mass accretion rate in units of 
10 _8 MQ2/r _1 and f is the stellar rotation velocity ex- 
pressed as a fraction of the star's breakup velocity (~ 
.15 for the CI/FS objects reported here). The large mass 
accretion rates believed to exist during the CI/FS stage 
(M Q ~ IQ-^M^yr- 1 ; Greene & Lada 2002) thus pre- 
dicts a disk-locking timescale on the order of 10 4 years, 
an angular momemtum evolution timescale much shorter 
than that assumed in §5.2. This timescale is sufficient to 
quickly remove significant amounts of angular momen- 
tum from the protostar and also consistent with the de- 
crease in stellar rotation observed between the CI/FS 
and CTTS phases. The reduced mass accretion rates 
characteristic of the CTTS phase, however, suggest a 
corresponding increase in the disk-locking timescale and 
thus less efficient extraction of angular momentum, con- 
sis tent with the spin-up of the rotation rate observed 
by ISicilia-Aguilar et all l|2005fl between the CTTS and 
WTTS phases. A decline in the mass accretion rate of a 
protostar as it evolves through the Class I/CTTS/WTTS 
phases and the inverse dependence of the timescale of the 
disk locking phenomena on the mass accretion rate ap- 
pear to mesh well with current observations of pre-main 
sequence rotation. 

6. CONCLUSIONS 

We have analyzed the projected rotation velocities (v 
sin is) of 38 D05 Class I/flat spectrum sources in nearby 
star formation regions with no apparent close compan- 
ions, allowing a robust investigation of the angular mo- 
mentum content of heavily embedded protostars. The 
key results of this study are: 

1. Class I/flat spectrum objects rotate more quickly 
(median equatorial rotational velocities ~ 38 km 
s _1 ) than CTTSs (median equatorial rotational ve- 
locities ~ 18 km s _1 ). Viewed in context of the 
protostellar evolutionary paradigm whereby Class 
I/flat spectrum objects are the less evolved prede- 
cessors of the older CTTSs, it is hard to explain the 
change in rotation rates between the two phases 
without the presence of some mechanism which 
transfers angular momentum out of the protostar 
between the Class I and CTTS phases. 

2. Projected rotation velocity (v sin i) seems to be 
weakly correlated with T e f f in our sample. Pro- 
jected rotation velocity does not seem to correlate 
with Br 7 emission (a common accretion tracer), 
the amount of excess continuum veiling (r&), or the 
slope of the SED (a) between the near and mid IR 
(a). 

3. Assuming Class I/flat spectrum sources possess 
physical characteristics (M*,R*,B*) typical of pre- 
main sequence stars, fully disk locked Class I ob- 
jects should have co-rotation radii within their pro- 
tostellar disks that match well (within 30%) with 
the pr edicted magnetic coupling radii of iShu et al.l 
(I994J) for similar sources. 



4. The factor of two decrease in rotation rates between 
Class I/flat spectrum and CTTS sources, when in- 
terpreted in the context of disk locking models, im- 
plies a factor of 5 or greater drop in mass accretion 
rates between the two phases. 

5. A lower limit of M ~ 10 -8 M Q yr _1 is required 
for objects transitioning from Class I/flat spectrum 
objects to CTTSs to account for the difference in 
rotation rates of the two classes by angular mo- 
mentum extraction through a viscous disk and no 
changes in protostellar physical structure. This 
lower limit implies a rate properly intermediate to 
typical estimates of Class I accretion rates (10 -6 
Mq yr" 1 ) and CTTS accretion rates (10~ 8 M Q 
yr- 1 ). 



The authors wish to recognize and acknowledge the 
very significant cultural role and reverence that the sum- 
mit of Mauna Kea has always had within the indigenous 
Hawaiian community. We are most fortunate to have 
the opportunity to conduct observations from this moun- 
tain. All data have been reduced using IRAF; IRAF is 
distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
This research has made use of NASA's Astrophysics Data 
System Bibliographic Services, the SIMBAD database, 
operated at CDS, Strasbourg, Franc e, and the VizieR 
datab ase of astronomical catalogues l)Ochsenbein et alJ 
2000). K.R.C gratefully acknowledges the support of 
NASA grant 80-0273. All authors acknowledge support 
from NASA Origins of Solar Systems program via UPN 
344-39-00-09 while conducting this work. The authors 
thank the anonomyous referee for a quick, thorough and 
illuminating report whose comments have improved the 
work presented here. 



Angular Momentum of Protostars 



13 



REFERENCES 



Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, 
A&A, 337, 403 

Barnes, S. A., Sofia, S., Prosser, C. F., & Stauffcr, J. R. 1999, ApJ, 
516, 263 

Bodenheimer, P. 1995, ARA&A, 33, 199 

Bontomps, S., Andre, P., Kaas, A. A., Nordh, L., Olofsson, G., 
Huldtgren, M., Abergel, A., Blommaert, J., Boulanger, F., 
Burgdorf, M., Cesarsky, C. J., Cesarsky, D., Copet, E., Davics, 
J., Falgarone, E., Lagache, G., Montmerle, T., Perault, M., Persi, 
P., Prusti, T., Puget, J. L., & Sibille, F. 2001, A&A, 372, 173 

Bouvier, J. 1990, AJ, 99, 946 

Bouvier, J., Bertout, C, Benz, W., & Mayor, M. 1986, A&A, 165, 
110 

Briceno, C, Calvet, N., Kenyon, S., & Hartmann, L. 1999, AJ, 118, 
1354 

Briceno, C, Hartmann, L. W., Stauffer, J. R., Gagne, M., Stern, 

R. A., & Caillault, J. 1997, AJ, 113, 740 
Chandrasekhar, S. & Munch, G. 1950, ApJ, 111, 142 
Claret, A. & Gimenez, A. 1990, Ap&SS, 169, 215 
Clarke, C. J. & Bouvier, J. 2000, MNRAS, 319, 457 
Covey, K. R., Greene, T., Doppmann, G., & Lada, C. 2005, AJ, 

'in prep.' 

D'Antona, F. & Mazzitelli, I. 1994, ApJS, 90, 467 
de Jager, C. & Nieuwenhuijzen, H. 1987, A&A, 177, 217 
Doppmann, G. W., Greene, T. P., Lada, C, & Covey, K. 2005, AJ, 
submitted 

Doppmann, G. W., Jaffe, D. T., & White, R. J. 2003, AJ, 126, 
3043 

Durisen, R. H., Mejia, A. C, & Pickett, B. K. 2003, Recent 
Research Development in Applied Physics, 1, 173 

Edwards, S., Strom, S. E., Hartigan, P., Strom, K. M., Hillenbrand, 
L. A., Herbst, W., Attridge, J., Merrill, K. M., Probst, R., & 
Gatley, I. 1993, AJ, 106, 372 

Gaige, Y. 1993, A&A, 269, 267 

Gameiro, J. F. & Lago, M. T. V. T. 1993, MNRAS, 265, 359 

Gray, D. F. 1992, The observation and analysis of stellar 
photospheres (Cambridge Astrophysics Series, Cambridge: 
Cambridge University Press, 1992, 2nd ed., ISBN 0521403200.) 

Greene, T. P. & Lada, C. J. 1997, AJ, 114, 2157 

— . 2000, AJ, 120, 430 

— . 2002, AJ, 124, 2185 

Greene, T. P., Wilking, B. A., Andre, P., Young, E. T., & Lada, 

C. J. 1994, ApJ, 434, 614 
Haisch, K. E., Greene, T. P., Barsony, M., & Stahler, S. W. 2004, 

AJ, 127, 1747 
Hartmann, L. 2002, ApJ, 566, L29 

Hartmann, L., Hewett, R., Stahler, S., & Mathieu, R. D. 1986, 

ApJ, 309, 275 
Hartmann, L. & Kenyon, S. J. 1996, ARA&A, 34, 207 
Hartmann, L. & Stauffer, J. R. 1989, AJ, 97, 873 
Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377 
Herbst, W., Bailer-Jones, C. A. L., & Mundt, R. 2001, ApJ, 554, 

L197 

Herbst, W., Rhode, K. L., Hillenbrand, L. A., & Curran, G. 2000, 
AJ, 119, 261 

Jayawardhana, R., Mohanty, S., & Basri, G. 2003, ApJ, 592, 282 
Johns-Krull, C. M. & Gafford, A. D. 2002, ApJ, 573, 685 
Johns-Krull, C. M., Valenti, J. A., Hatzes, A. P., & Kanaan, A. 

1999, ApJ, 510, L41 
Jones, B. F., Fischer, D. A., & Stauffer, J. R. 1996, AJ, 112, 1562 
Kaas, A. A., Olofsson, G., Bontemps, S., Andre, P., Nordh, L., 
Huldtgren, M., Prusti, T., Persi, P., Delgado, A. J., Motte, 
F., Abergel, A., Boulanger, F., Burgdorf, M., Casali, M. M., 
Cesarsky, C. J., Davies, J., Falgarone, E., Montmerle, T., 
Perault, M., Puget, J. L., & Sibille, F. 2004, A&A, 421, 623 
Kenyon, S. J. & Hartmann, L. 1995, ApJS, 101, 117 
Kenyon, S. J., Hartmann, L. W., Strom, K. M., & Strom, S. E. 

1990, AJ, 99, 869 
Konigl, A. 1991, ApJ, 370, L39 
Luhman, K. L. & Rieke, G. H. 1999, ApJ, 525, 440 



Magazzu, A., Rebolo, R., & Pavlenko, I. V. 1992, ApJ, 392, 159 
Martm, E. L., Koresko, C. D., Kulkarni, S. R., Lane, B. F., & 

Wizinowich, P. L. 2000, ApJ, 529, L37 
Massey, P. 1997, A User's Guide to CCD Reductions with IRAF , 

National Optical Astronomy Observatory 
Massey, P., Valdes, F., & Barnes, J. 1992, A User's Guide to 

Reducing Slit Spectra with IRAF, National Optical Astronomy 

Observatory 

Mathieu, R. 2003, in Proc. of the IAU Symposium 215, A. Maeder 
and P. Eenens; Ed., in press (see astro-ph 03031991 

Matt, S. & Pudritz, R. E. 2004, ApJ, 607, L43 

McLean, I. S., Becklin, E. E., Bendiksen, O., Brims, G., Canfield, 
J., Figer, D. F., Graham, J. R., Hare, J., Lacayanga, F., Larkin, 
J. E., Larson, S. B., Levenson, N., Magnone, N., Teplitz, H., & 
Wong, W. 1998, in Proc. SPIE Vol. 3354, p. 566-578, Infrared 
Astronomical Instrumentation, Albert M. Fowler; Ed., 566-578 

Mohanty, S. & Basri, G. 2003, ApJ, 583, 451 

Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 2965 

Neuhaeuser, R., Sterzik, M. F., Schmitt, J. H. M. M., Wichmann, 
R., & Krautter, J. 1995, A&A, 297, 391 

Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 

Padgett, D. L. 1996, ApJ, 471, 847 

Popham, R. 1996, ApJ, 467, 749 

Preibisch, T. & Smith, M. D. 1997, A&A, 322, 825 

Press, W. H., Flannery, B. P., Teukolsky, S. A., & Vetterling, 
W. T. 1986, Numerical recipes: the art of scientific computing 
(Cambridge University Press) 

Prosser, C. F. 1992, AJ, 103, 488 

— . 1994, AJ, 107, 1422 

Queloz, D., Allain, S., Mermilliod, J.-C, Bouvier, J., & Mayor, M. 

1998, A&A, 335, 183 
Rebull, L. M. 2001, AJ, 121, 1676 

Rebull, L. M., Wolff, S. C, & Strom, S. E. 2004, AJ, 127, 1029 
Rhode, K. L., Herbst, W., & Mathieu, R. D. 2001, AJ, 122, 3258 
Shu, F., Najita, J., Ostriker, E., Wilkin, F., Ruden, S., & Lizano, 

S. 1994, ApJ, 429, 781 
Sicilia-Aguilar, A., Hartmann, L. W., Szentgyorgyi, A. H., 

Fabricant, D. G., Furesz, G., Roll, J., Conroy, M. A., Calvet, 

N., Tokarz, S., & Hernandez, J. 2005, AJ, 129, 363 
Simon, M., Ghez, A. M., Leinert, C, Cassar, L., Chen, W. P., 

Howell, R. R., Jameson, R. F., Matthews, K., Neugebauer, G., 

& Richichi, A. 1995, ApJ, 443, 625 
Sneden, C. A. 1973, Ph.D. Thesis 

Soderblom, D. R., Stauffer, J. R., Hudon, J. D., & Jones, B. F. 

1993, ApJS, 85, 315 
Stassun, K. G., Mathieu, R. D., Mazeh, T., & Vrba, F. J. 1999, 

AJ, 117, 2941 

Stassun, K. G., Mathieu, R. D., Vrba, F. J., Mazeh, T., & Henden, 

A. 2001, AJ, 121, 1003 
Stauffer, J. R., Hartmann, L., Soderblom, D. R., & Burnham, N. 

1984, ApJ, 280, 202 
Stauffer, J. R. & Hartmann, L. W. 1987, ApJ, 318, 337 
Stauffer, J. R., Hartmann, L. W., & Jones, B. F. 1989, ApJ, 346, 

160 

Stauffer, J. R., Hartmann, L. W., Prosser, C. F., Randich, S., 

Balachandran, S., Patten, B. M., Simon, T., & Giampapa, M. 

1997, ApJ, 479, 776 
Tout, C. A. & Pringle, J. E. 1992, MNRAS, 256, 269 
Walter, F. M., Brown, A., Mathieu, R. D., Myers, P. C, & Vrba, 

F. J. 1988, AJ, 96, 297 
White, R. & Hillenbrand, L. 2004, ArXiv Astrophysics e-prints 
White, R. J. & Basri, G. 2003, ApJ, 582, 1109 
White, R. J., Ghez, A. M., Reid, I. N., & Schultz, G. 1999, ApJ, 

520, 811 

Wichmann, R., Torres, G., Melo, C. H. F., Frink, S., Allain, S., 
Bouvier, J., Krautter, J., Covino, E., & Neuhauser, R. 2000, 
A&A, 359, 181 

Wilking, B. A., Lada, C. J., & Young, E. T. 1989, ApJ, 340, 823 



